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An algorithm is presented that combines the ME+PS approach to merge sequences of tree- 
level matrix elements into inclusive event samples pQ with the POWHEG method, which 
combines exact next-to- leading order matrix elements with parton showers [2], The quality 
of the approach and its implementation in Sherpa [3] are exemplified by results for e + e~ 
annihilation into hadrons at LEP, for Drell-Yan lepton-pair production at the Tevatron 
and for Higgs-boson and VK + VK~-production at LHC energies. 

1 Introduction 

Facing the huge progress at the LHC, with first data taken, and first results already published, 
it is crucial to have reliable tools at hand for the full simulation of Standard Model signal and 
background processes as well as for the simulation of signals for new physics. This task is 
universally handled by Monte-Carlo event generators like Sherpa [3]. 

One of the key features of such advanced Monte-Carlo programs is the possibility to consis- 
tently combine higher-order tree-level matrix element events with subsequent parton showers 
(ME+PS) [T|. This feature has proved invaluable in various recent analyses of data from 
previous experiments, which arc sensitive to large-multiplicity final states. Despite being a 
tremendous improvement over pure leading-order theory, ME+PS merging still suffers from 
one major drawback of all tree-level calculations, which is their instability with respect to scale 
variations. This deficiency ultimately necessitates the implementation of NLO virtual correc- 
tions in Monte-Carlo programs. Two universally applicable methods were suggested in the past, 
which can perform this task, and whereof one is the so-called POWHEG algorithm [2]. This 
technique has been reformulated in [4], such that it can be applied in an automated manner. 

Having implementations of both, ME+PS merging and the POWHEG method at our dis- 
posal, the question naturally arises, whether the two approaches can be combined into an even 
more powerful one, joining their respective strengths and eliminating their weaknesses. A first 
step into this direction was taken independently in [5] and in [BJ- Here we will summarise the 
essence of the algorithms presented ibidem and exemplify the quality of related Monte-Carlo 
predictions. 
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2 The MENLOPS approach 



A formalism allowing to describe both, the ME+PS and the POWHEG method on the same 
footing was introduced in 0]. To compare, and, ultimately, to combine both methods, only the 
expressions for the differential cross section describing the first emission off a given core process 
must be worked out; this is where the combination takes place. 

In a simplified form, the expectation value of an observable in the POWHEG method can 
be described by the following master formula (for details see OH]) 



(O) 



POW 



=J2 J^bM^b) 



A i (t )O($ B ) + X! / d ^R\B^4^ 



Bi($ B ) 



.no omission 



(1) 



Dived 



omission 



where Bi($s) is the NLO-weighted differential cross section for the Born phase-space configu- 
ration $ B and Ai(i) = exp j - £\ J t d$ R \ B Rj/Bj | is the so-called POWHEG-Sudakov form 
factor. The indices i and j label parton configurations, see [4]. The parameter t is the ordering 
variable of the underlying parton-shower model and to is the respective cutoff. Hence, t is one 
of the variables used to parametrise the radiative phase space $> r \b- 

In a similar manner, a simplified master formula for the expectation value of O in the 
ME+PS approach can be derived. It reads (for details see 015]) 
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The terms labelled "ME domain" and "PS domain" describe the probability of additional 
QCD radiation according to the real-radiation matrix elements and their corresponding parton- 
shower approximations, respectively. In this context, K, R .\B i are the parton-shower evolution 
kernels and C Rj / Ri are the parton luminosities of the real-emission and the underlying Born 
configurations. In contrast to Aj(t) in Eq. (Q}, Aj(i) is the uncorrected Sudakov form factor of 
the parton-shower model. 

Combining the ME+PS method with POWHEG essentially amounts to combining the two 
above equations into a new master formula for the MENLOPS approach. This expression reads 
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In order to restore the POWHEG master formula, the "ME domain" term would have to be 
multiplied by the ratio of Sudakov form factors Aj(i)/Aj(i) only. Expanding this ratio to 
first order reveals that the above formula automatically yields next-to-leading order accurate 
predictions for any infrared and collinear safe observable O. 
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3 Results 



In the following we present selected results obtained with an implementation of the MENLO- 
PS algorithm in the Sherpa event generator. In particular we aim at detailing the improved 
description of data collected in various collider experiments. 

We focus first on electron-positron annihilation into hadrons at LEP energies (y/s =91.25 
GeV). Virtual matrix elements were supplied by BlackHat [7]. Figure [1] displays distributions of 
selected angular correlations in 4-jet production, that have been important for tests of pertur- 
bativc QCD. The good fit to those data proves that correlations amongst the final-state partons 
are correctly implemented by the higher-order matrix elements in the MENLOPS method. 

Similar findings apply in the analysis of the Drell-Yan process at Tevatron energies (y/s =1.96 
TeV). Figure [2] shows the transverse momentum distribution of the reconstructed Z-boson and 
the multiplicity distribution of accompanying jets, constructed using the D0 improved legacy 
cone algorithm with a cone radius of R = 0.5, p± t j > 20 GeV and \r]j\ < 2.5. The agreement of 
the MENLOPS result with the respective data is outstanding. 

We finally present some predictions for the production of Higgs-bosons through gluon-gluon 
fusion and for the production of >e + v e ] W~ [— > [i~v^\ at nominal LHC energies (yfs =14 

TeV). Virtual matrix elements for these analyses have been taken from [IT] and [12], respectively. 
Results are shown in Figs. [3] and 0] We observe very small uncertainties related to the intrinsic 
parameters of the MENLOPS approach. A detailed discussion is found in [5J. 
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Figure 1: The Korner-Schierholz-Willrodt (KSW) angle (left) and the modified Nachtmann- 
Reiter (NR) angle (right) in four-jet events defined using the Durham algorithm with y cut = 
0.008. Results at the parton level are compared to data from the OPAL experiment [8]. 
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Figure 2: Left: The transverse momentum of the reconstructed Z boson in Drcll-Yan events at 
the Tevatron (y/s = 1.96 TeV). Results obtained with the MENLOPS approach are compared 
to data from the D0 experiment [9]. Right: Inclusive jet multiplicity in Drell-Yan events. 
Monte-Carlo predictions are compared to data from jlO) . 
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Higgs boson p± Separation between leading and second leading jet 
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Figure 3: Left: The transverse momentum of the reconstructed Higgs boson in the gluon fusion 
process at nominal LHC energies (14 TeV). Right: Separation in rj-cj) space of the first and 
second hardest jet in Higgs-boson production via gluon fusion at nominal LHC energies. 
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Figure 4: Left: The transverse momentum of the reconstructed W + W~ system in > 
e + v e ] W~{-^ [i^Dfj,] events at nominal LHC energies (14 TeV). Right: Separation in ?y-0 space 
of the first and second hardest jet in W + W~ production at nominal LHC energies. 
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